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Abstract 

In general relativity, only relative acceleration has an observer-independend meaning: curva- 
ture and non-gravitational forces determine the rate at which world lines of test bodies diverge 
or converge. We derive the equations governing both in the conventional geometric formalism 
as well as using the background field method. This allows us to generalize the results to test 
bodies with charge and/or spin. The application of the equations to the motion of particles 
in a central field results in an elegant, fully relativistic version of the Ptolemaean epicycle 
scheme. 



1. World line deviation equations 

According to the equivalence principle, structureless test bodies (sometimes 
referred to as point masses) in a gravitational field move on geodesies of space- 
time. Their worldline x^(t) is a solution of the geodesic equation 

D 2 x^ _ d?x^ tl dx x dx u _ 

~ ~dP> + A - ~dT ~dT ~ °' (1) 

where the world-line parameter r is to be taken as proper time. Introducing the 
four- velocity as the time-like tangent unit vector to the worldline: = dx^/dr, 
the equation can be written in geometrical language as 

u-Vn = 0, u 2 = -1. (2) 

with V the covariant derivative. It is easily observed from eq. (|l|), that the 
proper acceleration = dPx^/dr 2 is not a covariant object. In particular, its 
vanishing or non- vanishing has no observer-independent meaning. 

In contrast, the relative acceleration between worldlines is a covariant quan- 
tity, and its vanishing or non-vanishing does not depend on the frame of ref- 
erence We recall the argument. Consider a one-parameter congruence of 
geodesies x M (r; A), where A labels the geodesies and r is the proper-time param- 
eter along the geodesic. We suppose the parametrization to be smooth, hence 
we can construct the tangent vector fields = dx^/dr, and = dx^/dX. It 
is straightfowardly established that 

As a corrolary, we obtain 

u-V(u-Vn) = u • V(n • Vu) = (u ■ Vn) • Vu + u^n u (V^V v u) 

= (n-Vu) ■Vu + u> i n v (V lt V v u) (4) 

= n • V(u • Vu) + u^n u [V^, V u ]u = u^n v 'R^[u, ■]. 
In component notation this reads 

D 2 ^ 



R K x/u K u»n\ (5) 



Dt 2 



The interest in the deviation vector obviously derives from the fact that, if 
Xq (r) = x^(t; Ao) is a solution of the geodesic equation (jl]), then to first order 
x i = x o + n^AX is a solution as well: 

3x^ 

x»(t; AO = x^(t; A ) + AA — (r, A ) « a^(r; A + AA). (6) 

oX 

It follows, that eq.(||) describes the covariant relative acceleration between these 
world lines. Of course, is only a first approximation to the neighboring 
geodesic at Ai = Aq + AA. To increase the precision of the approximation, 
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one has to compute higher-order derivatives w.r.t. A, by solving higher-order 
versions of eq.(||), involving not only the Riemann curvature tensor, but its 
derivatives as well. A systematic procedure of this type has been developed in 
ref . H . Here I pursue the first-order equation @ and study some generalizations 
and applications. 

We first observe, that eq.(|J) is linear and homogeneous in n^. It is therefore 
not very difficult to construct an action from which it can be derived. The 
lagrangean of interest reads 

r . . 1 Dn^ Dn v 1 „ A u „ 

In this lagrangean the metric, connection and curvature are those on the given 
reference geodesic Xq(t), with u^{t) = Xq representing the four- velocity along 
this same geodesic. These quantities act as background variables. Only the 
n^(r) are independent lagrangean generalized coordinates which are to be varied 
in the action. 

The action (@) can be derived independently by starting from the geodesic 
lagrangean 

and expanding x^{t) near the given background geodesic solution in the form 
x^ = Xq + n^AA. The term independent of A A does not contain and 
contributes a constant to the action. Next all terms linear in AA drop out 
of the result because xq is a solution of the geodesic equation. Finally, the 
terms quadratic in AA reproduce the expression (Q), up to a total proper- 
time derivative and terms which vanish because of the geodesic equation for 
Xq(t). Thus the lagrangean ([?]) represents the lowest-order non-trivial term in 
a systematic expansion: 

S[x] =mj cItL(x ) + ?n(AA) 2 J drL(n) + C[(AA) 3 ] (9) 

Obviously, the higher-order approximations can also be derived in this way. 

Clearly, as the variation of L(n) w.r.t. reproduces the geodesic deviation 
equation (g), the background- field method provides an alternative derivation of 
this equation. 

The above procedures can be generalized quite straightforwardly to cases in 
which test bodies are not completely structureless point masses, but carry e.g. 
charge and/or spin. In these cases particles do not move on geodesies, but on 
more general world lines ||, |J. For the case of charged particles in a combined 
electro-magnetic and gravitational field, the resulting world line deviation equa- 
tion was derived along the lines of eqs.@-@ in ref.[||]. An alternative derivation 
using the background field method starts from the action 

S g [x] = Jdr (^g^(x)x^ + qA^x)xA , (10) 
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with the overdot the usual short-hand for proper-time derivatives. The world- 
lines given by the stationary points of this action are solutions of the Einstein- 
Lorentz equation 

Dt 2 m u dr " 1 ' 

Now given a solution Xq(t) of this equation, and expanding the path in S q [x] 
as 

a^( T ) =^(r) + AA n^(r), (12) 
the action can be expanded to second order in AA as 

or! or! ( AA ) 2 /* j [" / Dn^ Dn u _ K A u „\ 

S q [x] = S q [x ] + ^^- J dr mlg^ — — + R^ KvX u K u x n»n v j 



+ 0[(AA) d ]. 

(13) 

To this order we then find that other solutions of the world- line equation ([Tl|), 
close to Xq(t), are given by (^), with the solution of the world- line deviation 
equation ||, ||] 

= Rxv^n" + — F* v + — V A n u"n\ (14) 

The alternative interpretation of n^, as parametrizing the distance between two 
particles on neighboring world lines, holds in this well, provided the par- 

ticles have the same charge-to-mass ratio q/m. Interestingly, this observation is 
not contradicted by the fact that one can obtain the Einstein-Lorentz equation 
as well as the world-line deviation equation ( |l4|) from reduction of the geodesic 
equation and geodesic deviation equation in five-dimensional space-time, as 
particles with different charge-to-mass ratio in four dimensions correspond to 
particles with different momentum in five-dimensional space-time ||. 

Similarly, pseudo-classical spinning particles can be described by the super- 
symmetric lagrangean ||, ||] 

1 i Dtp a 

L spin (x,f) = -g liu x^x +2^ a_ D7 r ' v 15 ) 

with an anti-commuting tangent-space vectorp] such that the pseudo-classical 
spin is described by S ab = — iift°"ip . The corresponding equations of motion for 
spinning particles can be written as 

D 2x>1 _ 1 cafe dm * DSab _ n / lfi \ 

Starting from a one-parameter congruence of solutions (a;^(r; A), ip a (r; A)) we 
define the deviation vectors 

d\> * dx ox » bV ' [ ' 



lr The transition between base-space and tangent-space vectors is made as usual by the 
vierbein e„ a and its inverse. 
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where uj^ a b is the spin connection. The covariant change in the spin-tensor is 
then 

n C*ab 

r b = _±- = -i + eV) ■ (is) 

These vectors satisfy the world-line deviation equations 



+ ~ (S a6 V A i?% u<V + r b R» uab if) . (19) 



DJ 



ab 



Dt 

They define the stationary points of the quadratic deviation action 

r I *s 1 Dnfl DnU i ^ Di a 1 k, \ a v 

L S pin{n,£) = 2 9 > J - I/ ~7)~~£) h 2 ^ a ~Dr + 2 ^ KvXU unn 

- % - ^ U^ab n "^ + V^xab u x n^ ( 20 ) 

2. Application: the Coulomb-Reissner-Nordstrom field 

World-line deviation equations can be used to compute the relative motion 
between particles in given background fields, or to obtain an approximation 
to solutions for orbits close to a known one. We illustrate the general results 
with an application to the study of the motion of charged particles in a central 
gravitational and electric Coulomb-Reissner-Nordstrom field. 

The vector potential and electric field strength for the Coulomb part of this 
solution of the Einstein-Maxwell equations are given by the one- and two-forms 

A = --Q-dt, F = dA = -Q-^dr Adt, (21) 
Amr 4-Kr 2 

whilst the metric for the gravitational field can be taken as 

- dr 2 = -B(r)dt 2 + — T dr 2 + r 2 (d6 2 + sin 2 6 dp 2 ) , (22) 
B(r) V / 

where B{r) = 1 — (2M/r) + (Q 2 /r 2 ); Q and M are the charge and mass of the 
central body which is the source of the field. 

The orbits of particles with mass m and charge q in this background can be 
computed in closed form in terms of elliptic integrals More precisely, the 
orbits are given by 

r ^ = n r ° < v ( 23 ) 

1 + e cos y((p) 

where y(<p) is the solution of the differential equation 

= JA + Bcosy + Ccos 2 y, (24) 
dcp v 
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with coefficients of given by 



Q 2 I ( q V\ 6M Q 2 . 9N 



£ 2 \ V 47rm / / ro 



o 



B = _~(m-^!], (25) 




C 



Here £ is the constant angular momentum per unit of mass. As the periastra 
of the orbit are at the points y(<p) = 2im, one can now compute the angular 
distance Ay between successive periastra. Writing Ay = 2tt + Sep, it follows 
that the periastron shift per orbit is 

»-*(--£-)+■- < 26 > 

y r 2Mr J 
the dots denoting terms of higher order in e, M/tq or Q/vq. 



Eqs.(23) and (Q) describe a general orbit in the exterior region of the cen- 
tral body. However, they do not provide all information about the orbit. In 
particular, as the time coordinate has been eliminated from these equations, the 
solution does not tell us where in its orbit the test particle is at any moment. 
Such information can be relevant for some important applications, e.g. to com- 
pute estimates of the amount of electro-magnetic and gravitational radiation 
emitted by the system. The method of world-line deviations is useful to obtain 
parametrized expressions of orbits (r(t),ip(t)). 

As the reference orbit, the zeroth order approximation to the real orbit, we 
take a circular one with constant radial coordinate R. Constants of motion on 
all orbits are the angular momentum per unit of mass, £ = u>R 2 , with oj = tp 
the angular velocity, and the energy per unit of mass e, defined by 

dt e — qQ/ AirmR 



dr 1-2M/R +Q 2 /R 2 ' 
On circular orbits the constants R, £ and e are then related by 



(27) 
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«-^b)' = 1 1 - : 5 l+ ^J1 1+ ^J- (28) 



and 




3f\ Q2/ 2_e_ N 

r 2 r r\ + r 2 




. (29) 

As all orbits are planar, we can always choose the orientation of the coordinate 
system such that 9 = n/2 for the reference orbit. For orbits tilted w.r.t. this 



one, we then find from eq.(14) that 

n e + to 2 n 9 = 0, (30) 
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from which it follows, as the physics dictates, that the distance perpendicular 
to the plane of the reference orbit oscillates with the period of the circular orbit 
T = 2-k/uj = 2ttR 2 /£. 

Considering orbits in the plane of the reference orbit, the world-line devia- 
tion equations (|i~4|) for the other components n l = (n t ,n r ,n ip ) become 



n { + 7 i ? n i + m^vP =0, (31) 
where the coefficient matrices take the form 



7 = l\ Y v \, rn = m r r . (32) 





This represents a system of coupled linear oscillators, which has solutions 

n*(r) = ngsinwiT, n r (r) = UqCosuiit, n^(r) = sinu^r, (33) 



where u\ is the solution of the characteristic equation for (pl|). The detailed 
form of this equation, using explicit expressions for the elements of the matrices 
7 and m were given in ||]. The resulting expression for the characteristic 
frequency is 

/ 3M Q 2 \ 

"> = "{ 1 --R + mR + ~)' (34) 

where the dots represent terms of higher order in M/R, Q/M or q/m. We also 
observe, that the amplitudes are not all independent: as u 2 = — 1 both on 
the original orbit and on the displaced world-line, it follows that n is space-like 
and u-n = 0. In the present case this general result translates to the constraint 

-4^)"»-4^" V »- £ "»= - < 35 > 

Although the components define the direction of the deviation, they do not 
determine the actual distance between neighboring world lines; this is given by 
equation (|i~2|) as Ax^ = n^AX. Therefore, for any particular orbit specified 
by the circular reference orbit (zeroth order approximation) and a world-line 
deviation vector n (first order approximation), we must determine in addition 
the scale factor AA to be applied. This can be done as follows. Comparing the 
approximate solution ( |3"I|) with the exact solution (^), we observe that 

r(<p) = R + Ar ^ R - eRcosy(ip). (36) 

Hence at the periastron, one has 

Ar = -AXn r = -eR. (37) 

Thus the scale is set by the eccentricity of the orbit. 
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Fig.l: Epicycle approximation to orbit with eccentricity e = \Ar\/R. 

Finally, we can determine the shift in angular coordinate between successive 
periastra, i.e. the advance of the periastron per orbital period. First observe, 
that the periastron occurs at the minima of ra r (r), i.e. for r n = (2n + 1)~k/oj\. 
Thus the amount of proper time elapsing between periastra is Ar = 2ir/u\\ the 
corresponding period of observer time is 

p2tt/u>i flf p2ir/ui 9-7r 

T= dr—= dr (u* + n*AA) = — u*. (38) 

Jo dr Jo uji 

Here u l is the rate of change of t per unit of proper time along the circular 
reference orbit. Next we observe, that at the proper times T n the angular 
coordinates at the reference orbit and the true orbit coincide: n if {r, n ) = 0. 
Hence the change in angular coordinate ip between successive periastra is the 
same as the change of this coordinate along the circular reference orbit after 
time T. This we can easily compute. Defining 

5<p = <p(to + T) — <p(t ) - 2vr, (39) 

and using for the angular velocity the expression d(p/dt = ipdr/dt = u/u t , we 
find 

ujT f uj \ ( 3M Q 2 \ 

Sip = — - 2tt = 2?r 1 « 2vr — . (40) 

Y u* \ut J \ R 2MR J v ' 

This is in perfect agreement with the expression (|26| ) obtained from the ana- 
lytical form of the orbit. 

It is of interest to consider the geometrical interpretation of the approxima- 
tion scheme we have used in a little more detail. The zeroth order approximation 
to the orbit we have constructed is a fully relativistic circular solution of the 
Einstein-Lorentz equation in a Coulomb-Reissner-Nordstrom field, with period 
To = 2tt/ujo- Included in this result is of course the simpler case of a circu- 
lar geodesic in a Schwarzschild field. The first-order correction is a geodesic 
deviation which oscillates in all its components in the same plane with period 
T\ = 2tt/u}\. Geometrically this represents another circular movement on the 
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background of the zeroth-order solution, i.e. an epicycle, with period slightly 
different from the zeroth-order approximation. This has two immediate con- 
sequences: the orbit becomes eccentric, and the period between extrema of 
the orbit differs from the period of the average (zeroth order) circular motion. 
This is in contrast with Newtonian gravity, where the periods are equal. Thus 
the extrema of the orbit (periastron and apastron) are shifted compared to the 
Newtonian approximation, by the amount predicted by the analytic description 
of the orbit. 

It can easily be shown Q], that higher-order world- line deviations all sat- 
isfy linear harmonic-oscillator type equations. Thus, computing higher-order 
corrections to our result amounts to the construction of higher-order epicycles. 
For the case of orbits in a central field, the method of world-line deviations then 
becomes a fully relativistic version of the Ptolemaean scheme |J , which differs 
genuinely from the standard post-Newtonian approximation scheme because it 
uses the eccentricity of the orbit and the quantities M/R, Q/M as expansion 
parameters, rather than v/c. As such this scheme offers an alternative to post- 
newtonian calculations of binary systems in a different physical regime, e.g. in 
the calculation of radiative effects. 

Acknowledgement 

I am indebted to Richard Kerner and Roberto Colistete for enlightening ideas 
and enjoyable cooperation. 



References 

[1] J.L. Synge, Relativity, the General Theory (Amsterdam, North-Holland; 
1960) 

[2] R. Kerner, J.W. van Holten and R. Colistete Jr., Relativistic Epicycles, 
Class. Quantum Grav. 18 (2001), 4725; e-print |a7Xlv:gr-qc/0102099| 

[3] J.W. van Holten, Proc. Sem. Math. Structures in Field Theories 1986-87, 
eds. E.A. de Kerf and H.G.J. Pijls (CWI Syllabus Vol. 26), 109 

[4] LB. Kriphlovich, Sov. Phys. JETP 69 (1989), 217 

[5] A. Balakin, J.W. van Holten and R. Kerner, Class. Quantum Grav. 17 

(2000) , 5009 

[6] R. Kerner, J. Martin, S. Mignemi and J.W. van Holten, Phys. Rev. D63 

(2001) , 27502 

[7] A. Barducci, R. Casalbuoni and L. Lusanna, Nuovo Cim. 35A (1976), 377 
[8] L. Brink, P. Di Vecchia and P. Howe, Nucl. Phys. B118 (1977), 76 
[9] Ptolemaios, Almagest (The mathematical syntax, ca. 145 AD) 



8 



